Vortex lattices in strong type-II superconducting two-dimensional strips 
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We show how to calculate semi- analytically the dense vortex state in strong type-II superconducting 
nanostructures. For the specific case of a strip, we find vortex lattice solutions which also incor- 
porate surface superconductivity. We calculate the energy cost to displace individual vortex rows 
parallel to the surfaces and find that this energy oscillates with the magnetic field. Remarkably, 
we also find that, at a critical field H* below H C 2, this "shear" energy becomes strictly zero for the 
surface rows due to an unexpected mismatch with the bulk lattice. 

PACS numbers: 74, 74.60.Ec, 74.76.-w 
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Despite increased technological interest in the trans- 
port properties of the vortex lattice state (VLS) [jjj in 
type-II superconducting nanostructures, the theoretical 
understanding of the interplay between the VLS and 
the interfaces is still mostly qualitative. While exper- 
iments are routinely done in nanostructures with vari- 
ous interface geometries ||-^] and explore all values of 
magnetic field H and temperature T, the order parame- 
ter, \l/(r), has been almost exclusively calculated in only 
two regimes, (i) For fields above the upper critical field 
H C 2, which separates the VLS from the metallic phase, 
solutions of the Ginzburg-Landau differential equations 
revealed long ago the existence of a highly localized or- 
der parameter at the surfaces of the sample which can 
survive up to a higher critical field H c3 = 1.695H c2 ||. 
Recently, experiments in novel engineered nanostructures 
Pj have spurred the interest in this surface superconduc- 
tivity and new calculations of upper critical fields have 
been done for different systems ||||. (ii) On the other 
extreme, ^(r) has been calculated for values of H just 
above the lower critical field H&, which separates the 
Meissner state from the VLS. The order parameter takes 
the form of a dilute VLS in which the inter-vortex dis- 
tance, ao, is large compared to the coherence length £ 
(typical size of the vortex cores). The properties of this 
dilute VLS in the presence of a surface parallel to H 
were analyzed in the past and numerous studies in 
thin films have been reported over the years j^j . 

The dilute VLS, which is correctly described in the 
London approximation where the vortex cores are ig- 
nored, only exists in a narrow range of H above H c \ for 
strong type-II (n 3> 1) superconductors Instead, 
a dense vortex state occupies most of the H — T phase 
diagram down to H w 0.3H C 2- In this regime, the vor- 
tex cores fill most of the space, interact strongly among 
themselves, and, most importantly, interact directly with 
the interfaces. To the best of our knowledge, only numer- 
ical simulations for nanostructures such as thin films or 



to H) |Tl|] have been reported in this regime. Our goal 
is to show that, for systems with simple geometries and 
simple orientations with respect to H, one can minimize 
semi- analytically the Ginzburg-Landau energy functional 
and obtain the full solution of the order parameter in the 
dense vortex regime. We illustrate the procedure for a 
disorder-free two-dimensional strip perpendicular to H. 
The order parameter comprises, generically, both a dense 
VLS and enhanced superconductivity at the surfaces. In 
addition, this procedure allows us to calculate the energy 
cost to displace individual rows parallel to the surfaces. 
This "shear" energy, E s , which is directly related to the 
measurable flow stress S, is, however, elusive in numer- 



ical simulations 1 1 1 . 

At the mean-field level (excluding thermal fluctuations 
B ) the order parameter can be obtained from minimiza- 
tion of the Ginzburg-Landau functional 
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slabs (placed parallel to H) |10| and strips (perpendicular 
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where G and G n are the Gibbs free energies of the super- 
conducting and metallic states, respectively. The term 
[— iKV — e*A(r)/c] 2 /2m* is the kinetic energy operator 
for Cooper pairs of charge e* = 2e and mass m* = 2m 
in a vector potential A(r) associated with the magnetic 
induction h(r). The parameters a and /3 have the usual 
meaning p2| . 

If, for simplicity, we restrict ourselves to two dimen- 
sions, k becomes effectively infinite and h(r) = H [ p^[ . 
In order to be able to find analytically the solutions which 
minimize Eq. ^ we first expand the order parameter as 
*(r) = J2 P , n C P ,™*p,n(r), where %, n (r) are the normal- 
ized eigenfunctions of the kinetic energy operator. The 
geometry of the interfaces determines the appropriate 
gauge choice and quantum number p (which is usually 
the linear or angular momentum); the microscopic details 
of the interfaces determine the boundary conditions for 
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the eigenfunctions . To obtain the order parame- 

ter in the dense vortex regime it is sufficient to consider 
an expansion of \&(r) in the lowest band (LB), n = [ p"4| |. 
The LB expansion allows us to write Eq. (Q) as 

N c 

G-G n = Y,u P \Cp\ 2 + 

p 

§ E c *p. °U c ps Cp, J d * K 

Pl,P2,P3,Pi 

where N c is the number of components, a p — a + e(p) 
is the condensation energy of the p component, and 
the second term represents the "interaction" between 
Cooper pairs. The interfaces are responsible for the non- 
uniformity of the condensation energy through the ki- 
netic energy, e(p), of the Cooper pair. — a p increases 
near interfaces with insulators, which favors surface su- 
perconductivity (see Fig. [l]), and shrinks near those with 
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FIG. 1. Cooper pair LB structure (in units of the cyclotron 
energy Tkj c /2) in a W = 16£ strip for different values of H 
(shifted for clarity). Notice the bending of the LB near the 
surfaces. Dashed lines represent —a. As H increases [from 
(a) to (e)] the number of components of the order parameter 
increases progressively. When —a lies below the LB in the 
center of the strip [panel (e)] only surface superconductivity 
survives. 

It is the p dependence of a p which, in principle, makes 
the minimization of Eq. (|J) non-trivial. In order to ac- 
complish this task one has to find: (i) the optimum N c , 
(ii) the quantum numbers {p}, and (iii) the complex coef- 
ficients C p . Some simple considerations (which are spe- 
cific for each geometry) allow us to solve this problem 
analytically. As an illustration we calculate the order 
parameter in an isolated strip. The dimensions of this 
ideal strip are L y — > oo, L z — > 0, and W, while H is 
chosen in the z direction. The appropriate LB eigen- 
functions are (in the Landau gauge) 3>fe(r) oc e lky x k (x), 
where k is the wavevector in the y direction and Xk(x) are 
nodeless functions which may be calculated numerically 



subjected to the boundary conditions of zero current on 
both surfaces. The integrals in Eq. || become 

I(ki,k 2 ,k 3 ,k 4: ) = (LyL z y 1 S kl +k 2 ,k 3 +k4 

J dxxk 1 (x)Xk 2 (x)xk 3 (x)XkAx) (3) 

and may also be calculated numerically. 

To find the lowest energy configuration we minimize 
(2jlhe Gibbs free energy in Eq. (|J) for each possible iV c 
separately (we will skip the trivial case iV c = 1, which is 
only relevant for very narrow strips with W ss £ fl^l): 

One vortex row (N c = 2): From the symmetry of the 
system one can consider the two components in a pair 
(k, —k) with k > [see, e.g., Fig. |l|(e)]. The symmetry 
also tells us that we can choose Ck = C- k - (Since the 
overall phase is irrelevant and the relative one only de- 
termines the position of the row in the y direction, we 
consider both coefficients to be real and positive.) There 
are only two types of interaction terms: (i) I(k, k, k, k) — 
I(—k,—k,—k,—k), which represent the interaction be- 
tween Cooper pairs occupying the same state [in short- 
hand, I(k)], and (ii) I(k, —k, —k, k) = I(k, —k, k, —k) = 
I(—k, k, —k, k) = I(—k, k, k, —k), which represent the in- 
teraction between Cooper pairs occupying different states 
[in shorthand, I(k, —k)]. The Gibbs free energy becomes 

G-G n = 2a k C 2 k +(3A k Ct, (4) 

where Ak — I(k) + 2I(k,—k). Straightforwardly, one 
obtains the minimal values of the coefficients, Ck — 
\J —ak/PAk, and the minimum Gibbs free energy, 

G — G n = —oq./(3Ak, for any pair (fc, — k). Knowing 
the band structure and the interaction integrals, the 
minimum-energy pair, (k, — k), can now be found, and 
\l/(r) takes the form of one centered vortex row as long 
as the two components overlap significantly For 
If > i; and H > H C 2 the order parameter is always 
formed by two components which do not overlap and 
constitute the usual surface superconductivity on both 
sides of the strip [see Fig. |l|(e)]. 

Two vortex rows (N c = 3): In addition to the pair 
(k, —ft) with positive, real coefficients Ck — C-k a third 
component can be considered at ft = [see, e.g., Fig. 
|l|(a)] with a complex coefficient |Co|e^°. The Gibbs free 
energy becomes 

G - G n = 2a k C 2 k + a \C \ 2 + ^{2A k C* + /(0)|C | 4 + 

4[2J(fc, 0) + cos (20 o )/(fc, 0, -k)]C 2 k \C \ 2 }, (5) 

where a new type of interaction terms, I(k, —k, 0,0) = 
/(-ft, ft, 0,0) = 1(0,0, ft, -ft) = 1(0,0, -ft, ft), appear 
[l(k,0, —k) in shorthand]. The last term in Eq. || repre- 
sents the "correlation" energy between rows and, regard- 
less of the magnitude of the coefficients, the minimum 
energy is obtained for 0o = t/2, i.e., when the two rows 
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are offset by half a row period, Oq/2. Upon minimization 
one obtains 



C k = 



l<%| = 



G — G„, — 



' agSk - a k I(0) 
(3[I(0)A k - 2S%] 

I 2a k S k - a A k 



I p[I(0)Ak - 2Sf] 

-2a 2 k I{&) - a\A k + 4a k a S k 
2/3[I(0)A k - 2S%] 



(0) 

(7) 
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for any trio (-fc, 0, fc), where S k = 2I(k, 0) - I(k, 0, -k). 

Three vortex rows (N c = 4) : The four components 
may correspond to pairs (k 2 ,— k 2 ) and (fci,— fci) with 
k 2 > fci > [see Fig. 0(b)]. We assume Cfc 2 = C_fc 2 
and ICfcJ = IC-feJ. Once again two of the phases have 
been set to zero (4> k2 and 4>- k2 ). The other two, 4> kl and 
4>~ kll must be chosen to minimize the energy. To make 
things analytically accessible we drop, for the moment, 
the weak correlation terms between vortex rows in the 
resulting Gibbs free energy: 

G - G n = 2a kl \C kl | 2 + 2a k2 \C k2 \ 2 + 

P(A kl \C kl | 4 + A k2 \C k2 1 + 4:S klyk2 \C kl \ 2 \C k2 1 2 ), 



where S kl ,k 2 
obtain 

\C kl \ = 

c k2 = 

G — G n = 



I(ki,k2) + I(ki,— fe)- Minimizing we 



la kl A k2 — 2a k2 S kl ^ k2 



A kl A k2 ] 



<ot k2 A kl — 2a kl S kl>k2 



V ^S 2 klM -A kl A k2 ] 

a kl Ak 2 + ct\ 2 A kl - 4:a kl a k2 S kuk2 
mS 2 klM ~A kl A k2 ] 



(10) 

(11) 
(12) 



for any set of four components (— &2, — fci, fci, k%). The 
correlation terms between adjacent rows, which ulti- 
mately determine the relative position between them, can 
be recast in the form 



8 k2Ml 2(i\C kl | 3 C fc2 /(-fci, fci, k 2 ) x 
[cos (2<p kl - (f)- kl ) + cos (20_ fcl - <p kl )} . 



(13) 



Notice that for k 2 = 3&i, the three rows have the same 
period a^. The minimum contribution of this term to the 
energy corresponds to 4> kl — 4>-k t — ft, which sets the 
offset between adjacent rows to Oq/2. The total energy 
can now be approximately obtained by subtracting the 
"lock-in" energy 4f3\C kl | 3 C fc2 /(-fci, k x , k 2 ) from Eq. © 
p6[ . If fc 2 7^ 3fci the choice of phases is irrelevant since 
the period of adjacent rows is different and they cannot 
fall into place. The above choice of phases gives the worst 
relative position between non-adjacent rows, i.e., aligns 
non-adjacent vortices. However, the term that accounts 
for such correlation is proportional to J(— k 2 , —fci, fci, k 2 ) 



and can be safely dropped since the overlap of X-fe2 and 
Xki is negligible in the expected minimal solution. 

For more than three vortex rows one adds more com- 
ponents [see panels (c) and (d) in Fig. [lj and follows 
similar considerations as for N c — 4: (i) to neglect corre- 
lations between non-adjacent rows, (ii) to minimize with- 
out the correlation terms for adjacent rows, and (hi) to 
subtract the lock-in energy at the end. One can^thus ob- 
tain analytical expressions for the C k s and G — G n for 
any number of rows from which the optimum set of quan- 
tum numbers {fc} can be found. We will only discuss the 
results in the following. Figure shows the number of 
vortex rows for strip thicknesses W = 12£ and W = 16£ 
as a function of H. As expected, the number increases 
in a step-wise manner. Figure [l] shows several snapshots 
of e(fc) for the W = 16£ strip for some characteristic val- 
ues of H. As long as —a or "chemical potential" for the 
Cooper pairs (dashed lines) remains above the LB in the 
center of the strip [panels (a) to (d)], the number of vor- 
tex rows increases with H . When —a lies below the LB 
in the center, i.e., when H > H c2 [panel (e)], the Cooper 
pairs can only nucleate near the surfaces (surface super- 
conductivity). For the narrower strip [Fig. 0(a)] the 
(lumber of vortex rows that can be accommodated be- 
fore the quenching of bulk superconductivity is logically 
smaller than for the wider one [Fig. ||(b)]. (To the nu- 
merical accuracy of our calculations all the rows forming 
the lattice disappear at H — H c2 in the limit W — > oo.) 
Peaks in the measured magnetization of thin films have 
already been associated with these type of transitions in 
the number of vortex rows (^Jl(| • 




FIG. 2. Number of vortex rows and E s as a function of H 
for a strip with (a) W — 12£ and (b) W = 16£. Shaded areas 
correspond to values of H with lock-in between the surface 
rows and the bulk lattice. Non-shaded areas to values of H 
without lock-in. 

In contrast to the numerical simulations HJITjJilL our 
analytical minimization allows us to calculate the en- 
ergy cost to "shear" or displace individual vortex rows in 
the direction of the surfaces. This energy E s is samplc- 
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dependent and proportional to the shear modulus cqq of 
the vortex lattice. To calculate E s one compares the min- 
imum energy with that obtained by switching the phases 
away from their optimum values so that vortices in adja- 
cent rows become aligned with each other. Dashed lines 
in Fig. U correspond to E s for the surface rows. When- 
ever a new row is added to the strip, E s experiences a 
sudden jump. This is because the addition of a new row 
"squeezes" all the others against each other, increasing 
their shear energy. On qualitative arguments |3 17 only 
one of the two possible lattice orientations is expected 
in this system. Although near the transition points the 
lattice is fairly distorted with ajj 7^ ao (ao corresponds 
to the perfect triangular lattice), the expected orienta- 
tion is still apparent in our results. The experimental 
value of the flow stress, which is proportional to E s , has 
been reported in Refs. || although instead of having the 
VLS confined in strips they had it "trapped" in disorder- 
free grooves etched in a dirty film. Oscillations in the 
flow stress were reported for the whole range of H and 
they were attributed to dislocations that appear when 
the periodicities of the groove superlattice and the VLS 
did not match. This picture does not consider the pos- 
sible overall distortion of the lattice. From our calcu- 
lations we conclude that, since the whole lattice recon- 
structs abruptly at each row addition, dislocations need 
not be invoked to explain jumps in the flow stress. Note 
that the boundary conditions that mimic the groove in- 
terfaces are different from the ones used here, but, as long 
as the vortex-interface interaction is sufficiently strong, 
this basic conclusion holds. 




FIG. 3. Order parameter for a W — 16£ strip at 
H — 0.9H C 2- Note the difference in period between the sur- 
face vortex rows and the bulk lattice. 

Finally, the most remarkable result to emerge from our 
calculations is the fact that E s goes strictly to zero at a 
critical field H* « 0.73iJ c2 as W -> 00. This is due to 
the fact that the surface rows become incommensurable 
with the bulk lattice at this critical field. Shaded areas 
in Fig. U correspond to the situation where lock-in takes 
place between the two surface vortex rows and the bulk 
lattice and, consequently, there is always a finite value 



of E s for such rows. For H > H* up to H C 2 the con- 
densation energy near the surfaces always overcomes the 
lock-in energy. The surface rows prefer to have a shorter 
period than the bulk ones despite of the fact that lock-in 
is no longer possible. Figure [5] shows the order parameter 
for H near H C 2 where the mismatch in period between 
surface rows and bulk lattice is apparent. This mismatch 
might reflect in a premature quenching of the flow stress 
measured in deeply etched NbGe grooves in experimen- 
tal setups similar to those of Refs. ||. The deep etching 
would create the necessary spatial modulation of the con- 
densation energy which triggers this phenomenon. 
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